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Abstract 


A  high  turbulent  intensity  combustion  chamber  has  been  designed  in  order  to  operate  with 
compressible  (0.3  <  M  <  0.5)  and  preheated  (T=500K)  air-methane  combustion.  These  conditions 
will  allow  the  investigation  of  different  flame  regimes;  most  notably  the  proposed  ‘Thickened 
Flame’  regime.  Initial  design  and  flow  validation  has  been  completed  in  order  to  begin  further 
experimentation.  A  10  kHz  Time  Resolved  Particle  Image  Velocimetry  (TR-PIV)  system  and  a  3 
kHz  Planer  Laser  Induced  Fluorescence  (PLIF)  system  have  been  integrated  with  the  system  in 
order  to  diagnose  the  flow  field  and  the  flame  respectively.  The  exhaust  and  chamber  cooling 
subsystems  were  designed  to  comply  with  safety  regulations,  and  the  control  systems  were  set  up 
in  a  way  that  allows  automated  (LabVIEW)  and  user  controlled  sequencing.  This  work’s  main 
purpose  is  to  characterize  and  map  the  flow  properties  at  maximum  flow  conditions  (M  >  0.3)  in 
order  to  map  the  following;  flow  structures,  kinetic  energy  at  different  length  scales,  velocity 
fluctuations  and  turbulent  intensity.  PIV  measurements  are  carried  out  for  three  different 
volumetric  flow  rates  (V  =  25  scfm,  55  scfm,  and  115  scfm)  and  three  different  grid  geometries 
for  each  (BR=  67%,  61%,  and  54%).  Upon  finding  that  the  main  flow  structure  consisted  of  jets 
caused  by  the  grid  at  higher  flowrates,  experiments  were  done  with  a  grid  BR  =  67%  which  was 
thinner  (50%  of  original  thickness)  than  the  previously  used  grids.  The  thinner  grid  produced 
structures  which  showed  flow  breakup  not  seen  previously  due  to  the  jets,  however  the  grid  BR  = 
61%  showed  best  overall  flow  structure.  The  experimental  results  allowed  the  velocity  fluctuation 
and  length  scale  terms  to  be  quantified  and  can  be  used  with  future  flame  studies  to  determine 
flame  regime  location  according  to  the  Borghi-Peters  diagram. 
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Chapter  1:  Introduction 


The  fundamental  understanding  of  turbulent  premixed  combustion  is  one  of  the  major 
research  topics  of  interest  in  the  industry  and  it  has  become  one  of  the  prior  research  interests  of 
the  Air  Force  Office  of  Scientific  Research  (AFSOR).  Many  experiments  aim  to  successfully  study 
flames  at  highly  turbulent  conditions  in  order  to  gain  a  better  understanding  in  turbine  design  and 
emissions  control.  There  have  been  many  studies  which  aim  to  provide  useful  data  in  this  topic, 
but  there  are  still  many  uncertainties  associated  with  turbulent  combustion.  Designers  still  rely  on 
the  design-test-modify  method  which  is  not  only  time  consuming  but  also  expensive.  It  is  because 
of  this  that  the  industrial  community  is  searching  for  a  reliable  and  effective  numerical  model  in 
order  to  design  next  generation  combustors,  and  these  models  need  a  detailed  experimental 
analysis  of  the  fundamentals  of  turbulent  flame  structure  and  a  wide  range  of  experimental  data  in 
order  to  be  valid  design  methods. 

The  development  of  high  repetition  rate  laser  diagnostic  techniques  such  as  Particle  Image 
Velocimetry  (PIV)  and  Laser  Induced  Fluorescence  (LIF)  enable  the  study  of  turbulent  flow  and 
turbulent  combustion.  These  non-intrusive  techniques  provide  a  way  to  study  the  actual  high  speed 
flow  with  minimal  effect  to  the  flow  field.  A  safe  and  reliable  high  intensity  combustion  system 
has  been  developed  in  order  to: 

Study  the  characteristics  of  premixed  flame  structure  at  compressible  and  high 
turbulent  intensity  conditions. 

Experimentally  establish  the  limits  of  flame  regime  boundaries  under  these 
conditions. 

Generate  experimental  data  through  the  use  of  simultaneous  PIV  and  PLIF. 
Generate  experimental  data  to  quantify  turbulent  flow  characteristics. 

The  main  goal  of  this  study  is  to  complete  the  characterization  and  mapping  of  the  fluid 
flow  at  several  flow  rates,  including  maximum  operating  conditions  (Mach  >  0.3).  System 
validation  has  been  conducted  and  the  next  objective  is  to  increase  the  operating  conditions  in 
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order  to  have  a  full  understanding  of  the  flow  interactions  inside  the  combustion  chamber. 
Different  flow  rates  and  grid  geometries  are  used  in  order  to  have  a  wide  range  of  data  sets. 

1.1  Objective 

The  Center  for  Space  Exploration  Technology  Research  (cSETR)  at  the  University  of 
Texas  at  El  Paso  specializes  in  propulsion  and  combustion  research.  Of  the  many  instruments 
available  in  the  center,  is  a  10  kHz  Particle  Image  Velocimetry  (PIV)  and  a  3  kHz  Planer  Laser 
Induced  Fluorescence  (PLIF)  system.  The  facility  also  has  3  kHz  high  speed  imaging  system  and 
a  compressor  system  which  can  supply  190  scfm  of  air  at  a  supply  pressure  of  140  psi.  This  specific 
experiment  is  supported  by  the  Army  Research  Office  (ARO)  and  a  collaboration  with  Princeton 
University.  The  main  objectives  of  this  project  are  proposed  as  following: 

To  study  the  global  characteristics  of  a  backward  facing  step  stabilized  flame  in 
compressible  and  high  intensity  turbulent  flow  using  high  speed  flow  visualization 
technology. 

To  study  flow  and  flame  characteristics  using  Time  Resolved  Particle  Image 
Velocimetry  (TR-PIV)  and  OH/CH  Planer  Laser  Induced  Fluorescence  (PLIF) 
imaging. 

To  generate  components  of  turbulent  intensity  and  turbulent  kinetic  energy. 

To  define  the  effects  of  bulk  velocity  parameters,  turbulence  parameters  and 
recirculation  zone  characteristics  on  flame  location,  shape,  fluctuation,  reaction 
zone  and  local  displacement  speed. 

The  preliminary  tasks  of  designing  and  validating  the  system  have  been  completed,  and 
the  main  objectives  can  now  be  completed.  At  this  stage,  the  main  objective  is  to  do  a  complete 
study  of  the  flow  characteristics  using  the  PIV  system.  Various  flow  rates  and  grid  geometries  are 
studied  and  presented  in  order  to  determine  which  turbulator  (grid)  generates  the  best  flow 
structures. 
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Chapter  2:  Theory  and  Literature  Review 

2.1  Flow  over  a  Backward-Facing  Step 

The  backward  facing  step  is  a  well-studied  phenomenon  in  fluid  mechanics,  and  its  flow 
characteristics  is  what  will  be  used  in  order  to  stabilize  the  highly  turbulent  flame.  The  fluid 
behavior  due  to  a  backward  facing  step  geometry  has  been  characterized  over  the  years  and  this 
fluid  behavior  can  be  seen  in  Figure  2.1.  Early  studies  were  conducted  using  laminar  flow  in  order 
to  develop  a  method  for  resolving  the  issues  of  a  backward  facing  step  through  the  use  of 
computational  methods  [1].  Later,  more  thorough  studies  were  perfonned  by  Biswas  et  al  [3]  and 
Annaly  et  al  [4]  to  further  the  understanding  of  the  backward  facing  step  different,  more  turbulent 
conditions.  Extensive  studies  were  also  conducted  by  Yuan  and  So  [2],  which  included  more 
characterization  parameters  in  a  simple  flow  over  a  backward  facing  step.  Several  fluid 
characteristics  which  occur  in  this  geometry  can  now  be  modeled  more  accurately  and  for  higher 
turbulent  conditions. 
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Figure  2.1:  Fluid  behavior  over  a  backward-facing  step  [2] 


Flow  simulations  for  the  present  study  were  completed  and  validated  in  previous  works 


[5].  A  Large  Eddy  Simulation  (LES)  model  was  used  in  order  to  identify  the  flow  structures 


(recirculation,  reattachment  length,  and  length  scales)  and  characteristics  for  this  geometry. 


Further  discussion  regarding  initial  flow  validation  is  found  in  Chapter  3. 
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2.2  Turbulent  Flow 


Turbulence  in  fluid  mechanics  refers  to  flow  which  is  unpredictable,  unsteady  and  chaotic. 
Turbulence  is  also  used  in  combustion  applications  since  it  offers  better  mixing  of  reactants.  It  is 
important  to  understand  Turbulence  in  order  to  correctly  characterize  and  map  the  flow  field  in 
this  project. 

The  most  basic  way  of  determining  the  turbulence  of  a  certain  flow  is  through  the  use  of 
the  Reynolds  number.  The  term  derives  from  fluid  mechanics  pioneer  Osborne  Reynolds  [7].  The 
Reynolds  number  represents  a  dimensionless  parameter  which  considers  the  ratio  of  inertial  to 
viscous  forces  and  serves  to  determine  the  dominant  force  that  influences  the  fluid’s  behavior. 
Equation  2.1  shows  the  mathematical  representation  of  this  term. 


V 


(2.1) 


In  this  equation  U,  L,  and  v  represent  the  characteristic  velocity,  characteristic  length,  and 
fluid  kinematic  viscosity  respectively  [6].  Different  ranges  in  the  Reynolds  number  determine 
whether  the  flow  is  considered  laminar  or  turbulent.  Simply  speaking,  the  faster  the  velocity  of  a 
fluid,  then  the  more  dominant  the  inertial  forces  are  compared  to  the  ‘slowing’  caused  by  viscosity. 
Table  2.1  shows  how  laminar  and  turbulent  flows  are  characterized  for  different  conditions. 


Table  2.1:  Reynolds  number  characterization  for  different  flows 


Reynolds  number 

Laminar 

Re  < 2300 

Transitional 

2300  <  Re  <  4000 

Turbulent 

Re  > 4000 

While  a  higher  Reynolds  number  generally  indicates  higher  turbulence,  the  transition 
region  serves  to  accommodate  cases  where  material  roughness  or  obstructions  in  the  flow  may 
affect  flow  development. 
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2.2.1  General  Description  of  Turbulence 


The  characteristics  and  aspects  of  turbulent  flow  are  described  in  thorough  detail  by  the 
work  of  Pope  [9].  As  mentioned  in  the  previous  section,  turbulent  flow  exhibits  unpredictability 
and  chaos.  Figure  2.2  shows  a  plot  of  the  velocity  at  different  time  intervals  of  a  turbulent  flow. 


Figure  2.2:  Plot  showing  variation  in  velocity  measurement  for  a  turbulent  flow  [9] 

In  order  to  characterize  the  characteristics  of  turbulence  in  the  flow,  the  velocity  must  be 
resolved  into  two  components;  mean  and  fluctuation  as  shown  in  Equation  2.2. 

u(t)  =  u  +  u'(t)  (2.2) 

u'(t)=u(t)-u  (2.3) 

The  mean  velocity  u  can  be  calculated  by  taking  an  average  velocity  measurement  for  a  certain 
time  interval  which  is  much  larger  than  the  turbulent  time  scale  [10].  From  the  velocity  fluctuation, 
it  is  then  possible  to  calculate  the  root  mean  square  (r.m.s.)  value  of  this  turbulence  component. 

urms  ~  sdev  ( u'(t ))  =  jvar  (it'(t))  =  jit'(t)2  (2.4) 
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Turbulent  Intensity 

Turbulent  intensity  is  a  parameter  used  to  quantify  the  turbulence  of  a  flow  field  and  is 
generally  expressed  as  a  percentage. 

/  _  Ur^  (2.5) 

U 

In  an  ideal  flow  (i.e.  non  turbulent),  the  Urms  value  would  be  0,  therefore  the  turbulent 
intensity  is  also  0%.  The  turbulent  intensity  in  high  speed  flows  is  generally  over  10%.  This 
variable  is  one  of  the  crucial  variables  in  mapping  the  flow  field  for  the  present  work. 

Turbulent  Kinetic  Energy 

Similar  to  kinetic  energy  (KE),  the  turbulent  kinetic  energy  (TKE)  is  a  measure  of  energy 
per  unit  mass  due  to  the  velocity  in  a  flow.  Though  in  order  to  calculate  TKE,  as  opposed  to  KE, 

the  fluctuation  component  of  velocity  is  used  as  follows- 

TKE  —  i  (ify?  +  ufi  ^  (2.6) 

Again,  the  r.m.s.  value  of  the  velocity  fluctuation  is  used  to  calculate  this  variable  [9]. 

Turbulent  Length  Scale 

When  dealing  with  turbulent  flow,  the  characterization  of  energy  modes  and  length  scales 
is  crucial  to  mapping  flow  fields.  As  flow  becomes  turbulent,  viscous  effects  bring  about  the 
formation  of  energy  dissipating  eddies.  These  eddies  provide  a  way  to  determine  the  length  scales 
in  a  flow.  The  largest  scale  eddy  is  one  which  carries  the  largest  amount  of  energy  in  a  flow,  and 
this  energy  is  dissipated  into  smaller  eddies  until  energy  is  simply  dissipated  as  heat.  Richardson 
and  Kolmogorov  introduced  and  refined  the  theory  of  what  is  known  as  the  ‘energy  cascade’  and 
‘Kolmogorov  hypothesis’  which  explain  this  phenomenon  [9].  The  energy  cascade  can  be  summed 
up  in  the  following: 

Big  whorls  have  little  whorls, 

Which  feed  on  their  velocity; 

And  little  whorls  have  lesser  whorls, 

And  so  on  to  viscosity  (in  the  molecular  sense)  [11] 
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Kolmogorov  then  introduced  a  method  to  quantifying  these  parameters  to  be  universal  and 
applicable  to  all  flows.  The  energy  dissipation  starts  at  the  largest  scale  L,  to  the  smallest  scale  q 
as  shown  in  Figure  2.3.  Once  energy  reaches  the  dissipation  range,  the  smallest  scales  of  length, 
time,  and  velocity  can  be  calculated  based  on  the  viscosity  of  the  fluid.  These  small  scale 
parameters  are  defined  as  follows- 


V  =  (7)1 

(2.7) 

(2.8) 

Uv  —  (sv  )1/4 

(2.9) 

e  =  Ha! 

L 

(2.10) 

Here  s  represents  the  energy  dissipation  which  scales  according  to  the  characteristic 
velocity,  u0  and  length  scale.  It  can  be  seen  that  the  smallest  scales  are  defined  by  the  kinematic 
viscosity  v  of  the  fluid  and  the  energy  dissipation  s.  Using  Equation  2.8,  we  can  find  the  time 
interval  used  to  the  average  the  velocity  in  Equation  2.2  and  2.3  in  order  to  correctly  calculate  the 
mean  and  fluctuation  velocity  components. 


Universal  equilibrium  range  Energy-containing 

range 

I 

Dissipation  range  I  Inertial  subrange 

I 

- r - 1 - 

if  Idi 

Figure  2.3:  Visual  representation  of  the  energy  dissipation  in  a  turbulent  flow 

Grid  Induced  Turbulence 

The  use  of  a  wire  mesh,  usually  referred  to  as  a  grid,  can  be  used  to  induce  turbulence  in  a 
cross-flow.  This  can  cause  the  flow  to  become  more  turbulent,  and  thus  is  a  common  method  to 
induce  turbulence  in  a  flow.  The  grid  creates  homogenous  turbulence  which  increases  with  time 


I 

I 
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by  t  =  x/Uo,  and  decays  along  the  x-direction.  The  use  of  a  grid  shows  that  the  r.m.s.  value  of  the 
x-component  of  the  velocity  is  10%  greater  than  the  y-direction  component  [9],  The  flow  begins 
to  exhibit  transitional  turbulence  at  the  point  where  the  grid  is,  and  gradually  develops  into  a  fully 
developed  turbulent  flow.  As  the  flow  progresses  axially  the  flow  begins  to  decay,  starting  with 
the  smallest  scale  eddies.  A  visual  representation  of  grid  induced  turbulence  can  be  seen  in  Figure 
2.4.  In  the  present  study,  three  different  perforated  plates  are  used  in  order  to  increase  the 
turbulence  and  effectively  increase  the  velocity  fluctuation  and  turbulent  intensity. 


Stage  (ii) 

<tvs — '  /T- 


Stage  (iii) 

.C  ~rs/' — 


”5) 


Fully  developed  turbulence 

Transition  to 

fully  developed 

Small-scale 

turbulence 

eddies  decay  first 

Figure  2.4:  Visual  representation  of  grid  induced  turbulence  in  a  duct 

2.3  Optical  Diagnostics 


Laser  imaging  techniques  provide  a  non-intrusive  way  of  measuring  flow  parameters, 
combustion,  emissions  and  all  their  individual  components.  Optical  diagnostics  began  as  early  as 
the  1800’s,  but  were  popularized  in  the  early  20th  century  and  are  now  widely  used  in  turbulent 
combustion  research  [14].  There  are  a  wide  variety  of  imaging  techniques  and  systems,  and  in  the 
present  study  a  PIV  and  LIF  system  are  to  be  used.  These  systems  are  to  be  run  separately,  as  well 
as  simultaneously. 


2.3.1  Particle  Image  Velocimetry  (PIV) 

Particle  Image  Velocimetry  (PIV)  is  a  way  to  map  flow  characteristics  such  as  velocity 
vector  fields,  vorticity,  strain  and  divergence  of  the  flow.  In  general,  a  PIV  system  consists  of  the 
following  three  important  components;  seeding  particles  or  other  flow  tracer,  a  high  powered  laser, 
and  a  CCD  camera.  The  works  of  Adrain  [15]  explain  the  working  of  a  PIV  system  and  is  therefore 
a  useful  resource  to  help  understand  some  basic  concepts.  PIV  systems  are  usually  complemented 
with  software  which  generates  all  the  desired  parameters  based  on  particle  movement.  It  utilizes 
the  captured  images  from  the  camera  to  determine  the  velocity  components  by  measuring  the 
calibrated  displacement  of  a  single  particle  over  a  specified  time  step.  Images  are  recorded  and 

paired  in  order  to  have  a  ‘first’  and  ‘second’  image  which  the  displacement  can  be  calculated  from. 

Ax 

u  =  T7  (2.11) 


v  = 


Ay 

At 


(2.12) 


PIV  accuracy  depends  on  successful  synchronization  of  image  capture  and  laser  firing  as 
well  as  proper  calibration  and  image  capture  quality.  The  correct  use  of  particles  is  also  paramount 
to  obtaining  useful  data.  Particles  should  be  sufficient  enough  to  be  traced  by  the  camera  from  one 
image  frame  to  the  next,  but  also  sparse  enough  to  not  over-saturate  the  image.  A  relation  between 
particle  and  fluid  properties  can  be  used  to  have  adequate  sizing  and  volumetric  concentration  of 
particles.  This  can  be  seen  by- 


PpTcdpVr 

18n8p 


<  1 


(4.2) 


Where  dp,  vr,  Sp,  and  pp  are  the  particle  diameter,  fluid  viscosity,  average  distance  between 
particles,  and  particle  density  respectively. 
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Figure  2.5  shows  how  data  is  obtained  from  using  a  PIV  setup.  As  shown  in  the  figure,  the 
laser  is  fired  into  an  optical  component  which  changes  the  beam  from  its  original  Gaussian 
distributed  beam,  to  a  light  sheet.  This  sheet  then  illuminates  the  particles  in  the  flow,  which  are 
then  recorded  with  the  CCD  camera,  and  finally  post  processed  using  an  image  processing 
technique  or  software.  It  is  to  be  noted  however,  that  the  light  sheet  generated  is  not  completely  2- 
dimensional,  but  the  image  capture  is.  The  light  sheet  is  not  a  sheet,  but  a  volumetric  space  which 
can  be  reduced  by  using  the  right  optics,  and  distancing  all  equipment  appropriately.  It  is  because 
of  this  that  the  image  must  be  focused  in  the  mid  plane  inside  the  light  sheet  volume  where  the 
laser  intensity  is  highest. 

2.3.2  Planar  Laser  Induced  Fluorescence  (PLIF) 

The  works  of  Grisch  and  Orain  [17]  show  a  general  overview  of  laser  induced  fluorescence 
(LIF).  The  main  principle  behind  LIF  is  the  use  of  pulsed  laser  energy  to  increase  the  energy  in 
combustion  species  such  that  they  reach  what  is  referred  to  as  the  excited  state.  Once  the  molecule 
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reaches  this  excited  state,  then  it  will  begin  to  emit  radiation  energy  causing  it  to  fluoresce  until  it 
then  reaches  a  state  of  decay.  Figure  2.6  shows  a  clear  representation  of  what  happens  when  a 
particle  absorbs  enough  energy  from  the  laser. 


Figure  2.6:  Visual  representation  of  molecular  activity  after  energy  absorption  [17] 

This  allows  the  image  capture  of  the  molecule  using  a  CCD  camera  and  intensifier.  The 
fluorescence  is  captured  at  specific  wavelengths  according  to  different  species  as  seen  in  Table 
2.2.  These  wavelengths  provide  information  on  what  frequency  the  pulse  laser  should  be  fired  to 
cause  the  fluorescence  of  individual  molecules.  A  band  pass  filter  is  used  to  isolate  the  image 
capture  to  a  specific  wavelength  and  eliminate  any  imaging  caused  by  chemiluminescense.  The 
most  commonly  analyzed  combustion  radicals  are  OH,  CH,  and  CH2O.  These  molecules  can  be 
used  to  define  the  flame  front  in  different  ways.  OH  is  the  most  abundant  radical  present,  and  it 
also  offers  the  most  intense  image  capture  [18].  Most  studies  using  OH-PLIF  imaging  are 
successful,  and  this  molecule  is  heavily  favored  when  conducting  LIF  experiments.  One  limitation 
to  analyzing  OH  is  the  fact  that  this  radical  is  not  a  clear  indicator  of  the  exact  flame  front  position, 
or  flame  thickness.  This  is  due  to  the  fact  that  OH  can  generate  high  intensity  fluorescence  after 
the  flame  front.  CH2O  (formaldehyde)  is  also  be  used  to  generate  flame  imaging,  but  similarly  to 
OH,  CH2O  can  generate  fluorescence  outside  of  the  flame  front.  The  most  precise  method  for 
imaging  is  using  the  CH  radical.  This  radical  occurs  along  the  flame  front,  and  thus  can  provide 
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an  accurate  measure  of  the  flame  front  thickness  if  successful  imaging  is  obtained  as  shown  in 


Figure  2.7. 


Figure  2.7:  OH,  CH,  and  CH2O  imaging  for  a  simple  jet  [19] 


Dye  Lasers  are  used  due  to  their  ability  to  be  tuned  to  a  broad  range  of  wavelengths, 
particularly  in  the  visible  spectral  range,  the  near  ultraviolet  and  near  infrared  regions.  This  non- 
intrusive  optical  technique,  like  the  PIV,  resolves  important  information  about  the  structure  of  the 
flame  and  the  combustion  process  for  the  next  generation  combustor  design.  This  type  of  laser 
system  can  be  used  to  gather  experimental  combustion  data  which  can  be  used  to  validate  turbulent 
modeling. 


Table  2.2:  Absorption  wavelengths  for  different  combustion  species 


Molecule 

Absorption  wavelength  (nm) 

OH 

306 

C2 

516 

CH 

431 

CH2O 

320-360 

NO2 

450-470 

Appendix  II  shows  the  specifications  of  the  LIF  system  and  intensifier  used  for  this  proj  ect. 
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2.4  Flame  Regimes 


Flame  characteristics  change  in  accordance  to  the  fluid  conditions  during  combustion. 
There  are  various  methods  and  theories  explaining  these  characteristics,  but  a  widely  used  method 
for  classifying  these  different  types  of  flames,  was  proposed  by  Borghi,  and  later  modified  by 
Peters  [20-21]  ,  and  is  known  as  flame  regime  diagram  (Figure  2.8).  The  diagram  is  characterized 
by  dimensionless  parameters  (ReT,  Da,  Ka),  turbulent  parameters  (u1,  Lt),  and  flame  parameters 
(Sl,  5l).  Figure  2.9  shows  theoretical  pictorial  representations  of  these  flame  regimes,  and  how  the 
flame  fronts  are  affected. 
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Figure  2.8:  Flame  regime  diagram  by  Borghi  and  modified  by  Peter 


The  different  flame  regimes  are  explained  in  order  to  clarify  the  flame  front  behavior  at 
varying  turbulence.  The  region  of  interest  in  this  project  is  the  thin  reaction  zone  which  causes  a 
thickened  flame.  This  flame  regime  was  proposed  and  studied  by  Zimont  [22],  and  has  been 
studied  extensively  in  order  to  understand  turbulent  mechanics  in  combustion.  This  flame  regime 
can  be  seen  in  application  where  turbulence  is  a  determining  factor  in  design  such  as  in  turbines, 
ramjet,  and  scramjet  engines. 
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Figure  2.9:  Different  flame  regimes  and  the  graphical  illustrations  of  flame  front  characteristics 

Wrinkled  Flame  Regime 

The  wrinkled  flame  regime  occurs  when  flame  speed  is  higher  than  the  fluctuations  in  the 
flow  field,  and  thus  the  combustion  kinetics  have  a  higher  effect  on  the  flame  front.  The  flame 
thickness  is  smaller  than  the  integral  length  scale,  meaning  the  flow  eddies  are  much  larger  than 
the  flame  thickness  and  therefore  cause  a  turbulent  appearance  in  the  flame,  hence  the  wrinkled 
flame. 

Corrugated  Flame  Regime 

The  corrugated  flame  occurs  when  turbulent  intensity  predominates  over  flame  speed,  thus 
causing  breakup  in  the  flame  front.  Breakup  occurs  due  to  the  fact  that  Ka  <  1 ,  therefore  the  flame 
thickness  is  still  smaller  than  the  small  scales  in  the  turbulent  flow.  This  cause  the  flame  to  exhibit 
vortices  in  its  flame  front  which  can  cause  pockets  of  unburnt  reactants  to  be  present  in  the  flame. 

Thickened  Flame  Regime 

The  thickened  flame  occurs  when  the  turbulent  intensity  is  much  greater  than  the  flame 
speed,  causing  higher  levels  of  disruption  in  the  flame  front  than  in  the  previously  discussed 
corrugated  flame.  The  increase  in  turbulence  also  causes  the  small  scale  eddies  to  be  smaller  than 
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the  flame  thickness.  This  in  turn  causes  the  small  scale  eddies  to  enter  the  flame  reaction  zone  and 


cause  a  thickened  flame  front.  The  flame  is  also  stretched  and  the  reaction  zone  is  altered  when 
the  flame  is  thickened  such  that  the  flame  speed  can  increase  due  to  the  inner  reaction  zones  within 
the  flame  front.  This  present  study  system  aims  to  achieve  the  conditions  required  to  generate  a 
thickened  flame  experimentally.  In  this  work,  the  author  focuses  on  the  quantification  of  the 
turbulent  parameters  u’  and  Lt,  which  will  be  used  with  future  LIF  data  to  detennine  the  thickened 
flame  characteristics. 

Distributed  Flame  Regime 

This  flame  regime  occurs  when  the  chemical  kinetics  of  the  combustion  controls  the 
combustion  process.  A  distributed  flame  is  achieved  when  the  mixing  of  the  reactants  is  higher 
than  the  time  scales  during  combustion.  This  causes  a  flame  which  is  evenly  distributed,  and  whose 
flame  front  is  unlike  any  of  the  previously  mentioned  flames. 
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Chapter  3:  Experimental  Setup 


3.1  System  description 

The  design  of  this  turbulent  combustion  system  has  been  discussed  previously  [5].  As 
previously  mentioned,  the  current  system  is  to  be  operated  with  air  and  methane  as  reactants  in 
order  to  study  combustion  at  compressible  flow  rates  (M  >  .3).  The  system  requirements  which 
drove  the  design  of  this  system  can  be  seen  in  Table  3.1.  The  main  combustor  is  shown  in  Figure 
3.1,  and  consists  of  the  following;  combustion  chamber,  entrance  and  exhaust  regions,  mixing 
chamber,  quartz  windows,  and  modular  backward-facing  step.  The  other  subsystems  are  discussed 
in  detail  in  the  following  sections  to  provide  an  understanding  behind  the  overall  design  of  the 
high  turbulent  intensity  system. 


Table  3.1  System  Requirements 


System 

Requirement 

Combustion  chamber 

Optically  accessible 

Grid  turbulence  generator 

Changeable  step  size  (variable  dimension) 

Max  Pressure:  6  bar  (  87  psi) 

Reactants 

Air  /  CH4 

Air 

Pre-heated  at  400-600  K  (260  -  620  F) 

Velocity  Mach  >  0.3 

Flow  rate:  Max  at  0.4  kg/s 

CH4 

Flow  rate:  Variable  (max  0.07  kg/s) 

Inlet  pressures 

100  -  600  kPa  (14  -  87  psi) 

U'/Sl 

5-10 

Lt/5 

50-200 

An  exploded  model  of  this  subsystem  can  be  seen  in  Figure  3. 1,  and  a  cross-sectional  model 
showing  the  assembled  chamber  is  shown  in  Figure  3.2.  The  design  of  the  entrance  and  exhaust 
regions  is  based  on  a  polynomial  profile  curve  which  allows  for  the  entrance  and  exhaust  the 
combustor  to  follow  a  uniform  ‘top  hat’  velocity  profile  [23].  Extensive  FEA  and  CFD  modeling 
has  been  performed  on  this  subsystem  to  validate  the  design  and  ensure  it  can  meet  the  system 
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requirements,  while  also  maintaining  the  safety  of  those  operating  it.  The  FEA  analyses  were 
conducted  under  maximum  chamber  pressure  of  6  bar  and  maximum  temperature  of  500K.  Table 
3.2  shows  the  subsystem  components  and  their  respective  factor  of  safety  (FOS).  This  table  shows 
that  the  failure  mode  will  occur  at  the  windows;  more  specifically  the  top  window.  In  order  to 
address  this,  the  thermal  expansion  is  considered  when  choosing  a  gasket  for  the  combustor 
windows  and  flanged  unions.  The  chosen  gasket  material  is  a  homogenous  high  temperature 
application  graphite  gasket  sheet.  This  gasket  material  can  withstand  up  to  6000K  at  pressure  of 
10  bar,  which  is  well  above  the  maximum  conditions  the  system  will  operate  at.  This  gasket 
material  also  has  a  40%  compression  ratio  with  20%  recovery  factor,  which  offers  the  desired 
dimensions  once  thermal  expansion  is  considered  along  the  largest  stainless  steel  to  quartz 
interface.  A  perforated  plate  is  used  as  a  turbulator  in  order  to  cause  grid  induced  turbulence  as 
described  in  Chapter  2.  The  three  different  grids  consist  of  a  rectangular  plate  with  a  5x5  array  of 
holes.  The  different  grid  blockage  ratios  are  54%,  61%,  and  67%.  The  placing  of  this  grid  is  shown 
in  Figure  3.2. 


Mixing  Section 


Methane  Lines 


Air  Line 


facing  step 


Cooling 


Square  Inlet 


t 

Exhaust  Section 


Combustion  Chamber 


Figure  3.1:  Exploded  model  of  the  combustor  subsystem  and  its  components 
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Figure  3.2:  Cross-sectional  view  of  the  combustor  showing  2D  plane  used  in  simulations  [5] 


Table  3.2:  Components  in  combustor  and  their  factor  of  safety 


Component 

FOS 

Window  (thinnest) 

3.5 

Entrance/Exhaust  section 

11.8 

Grid 

16.8 

Combustor 

12 

The  main  air  line  consists  of  1”  316  Stainless  Steel  tubing  which  draws  air  from  a 
compressor  that  can  operate  at  190  scfim  and  140  psig.  This  flow  rate  allows  the  experiments  to  be 
conducted  at  compressible  air  conditions.  The  line  is  insulated  after  the  preheating  section  to 
prevent  temperature  drop  of  the  preheated  air.  Sizing,  pressure,  density,  and  Mach  number 
calculations  were  carried  out  using  compressible  flow  theory  [24],  An  overall  schematic  of  the 
system  can  be  seen  in  Figure  3.3.  Figure  3.4  and  Figure  3.5  show  the  pilot  flame  subsystem,  which 
will  serve  as  the  ignition  source  for  the  flame  in  the  combustor,  in  greater  detail.  This  pilot  flame 
subsystem  will  operate  using  two  hydrogen- air  premixed  flames  integrated  into  the  modular  step. 
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The  pilot  flame  is  ignited  using  performance  sparkplugs  with  modified  electrodes  which  cause  the 
electric  arcs  to  be  at  the  top  of  the  step  surface.  The  sparkplugs  are  powered  by  a  12V  lawn  mower 
battery  and  controlled  with  LS 1  coils  wired  to  a  TTL  5 V  signal  generator. 


Figure  3.3:  Flow  schematic  of  the  combustion  system  and  detailed  pilot  flame  flow  schematic 
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Figure  3.4:  Combustion  chamber  showing  location  of  pilot  flames 


Figure  3.5:  LSI  coil  wiring  diagram 
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3.2  Air  Preheating 


To  preheat  the  air  to  the  desired  maximum  500K  temperature,  a  large  heat  source  is  to  be 
considered.  The  constraints  for  determining  the  preheating  section  were  based  on  the  main 
laboratory  constraint  of  having  100A  of  current  with  a  3-phase,  480V  source.  The  most  effective 
way  to  utilize  the  electrical  power  to  heat  the  air  is  through  the  use  of  a  resistance  heater  with  high 
efficiency.  A  lOOkW  commercial  heater  with  a  10%  efficiency  will  be  used  to  preheat  the  air. 
Figure  3.6  shows  a  schematic  of  this  heater,  and  Figure  3.7  shows  the  heating  elements  encased 
within  the  heating  unit.  The  required  heat  vs  the  maximum  heat  that  can  be  obtained  was  first 
compared  in  order  to  ensure  the  possibility  of  preheating  the  air  to  the  required  temperature.  This 
was  done  via  a  simple  energy  calculation  required  for  the  air  and  maximum  electrical  power 
available. 

Q  —  m  *  cp  *  (Th  —  Tl)  (3.1) 

P  —  0.9  *  /  *  V  *  PF  *  cos (0)  (3.2) 

Where  the  air  properties  are  based  on  maximum  conditions,  V,  PF,  and  (j)  are  the  voltage, 
power  factor  (0.95  for  resistance  heaters),  and  phase  respectively.  These  calculations  show  the 
power  delivered  by  the  heater  is  sufficient  to  preheat  the  air  at  maximum  conditions. 


Figure  3.6:  Heater  (a)  model  with  stand  and  (b)  manufacturer  specific  drawing 
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Figure  3.7:  Inconel  heating  elements  inside  the  heater 
The  heating  elements  consist  of  36  hairpin  Inconel  elements  with  a  heat  flux  of  3.6  W/cm2. 
Using  this  information  and  assuming  convective  heat  transfer  between  the  air  and  heating 
elements,  we  can  find  a  value  for  the  exit  temperature  of  the  air. 

Talrlcut=Tin+^  (3.3) 

The  exit  temperature  of  the  air  is  5 15K,  which  gives  a  15  degree  margin  for  the  preheated 
air.  The  theoretical  value  was  then  validated  by  carrying  out  a  computational  model  using  the  full 
scale  dimensions,  material,  and  flow  properties.  A  transient  k-s  turbulent  model  was  used  in 
ANSYS  14.5  Academic  version  with  a  2D  mesh  consisting  of  a  400000  structured  elements  biased 
toward  the  exit  section  of  the  heater.  Viscous  heating  was  enabled,  and  second  order  upwind 
solutions  were  obtained.  The  transient  solution  is  carried  out  for  5  minutes,  then  10  minutes  in 
order  to  extract  a  steady  state  temperature  profile.  Figure  3.8  shows  the  time  resolved  solution, 
concentrated  at  the  exit  region  of  the  heater  and  its  temperature  contours.  Figure  3.9  shows  a  plot 
extracted  from  the  two  transient  solutions,  assuming  that  after  10  minutes  the  flow  is  considered 
steady  state.  Table  3.3  shows  a  comparison  between  the  theoretical  and  calculated  values  of  exit 
temperature.  It  should  also  be  noted  that  after  the  preheating  section  the  air,  though  insulated,  will 
lose  ~  10K  due  to  compressible  flow  effects.  Appendix  III  shows  the  compressible  flow  equations 
including  Fanno-line  flow  (assuming  isothermal  flow)  used  to  calculate  these  effects. 
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Figure  3.8:  Static  temperature  contours  of  preheated  air 


0  0.02  0.04  0.06  0.08  0.1  0.12 

Position  (m) 


Figure  3.9:  Plot  of  steady  state,  transient,  and  average  temperature  profiles  at  heater  exit 
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3.3  Exhaust  Design 


The  exhaust  design  was  developed  to  meet  the  temperature  requirement  of  exhaust  gases 
not  exceeding  temperatures  of  80  °C.  The  exhaust  fan  inside  the  laboratory  has  a  suction  capability 
of  230  scfm,  which  is  above  our  maximum  operating  conditions.  The  exhaust  temperature  of  the 
air  is  assumed  to  be  800K  at  maximum  operating  conditions  based  on  preliminary  combustion 
modeling  [5].  Several  methods  for  exhaust  cooling  were  proposed  including  the  use  of  a  coolant 
similar  to  the  combustor  cooling,  and  liquid  nitrogen  spraying  to  reduce  the  exhaust  gas 
temperature.  These  methods  however  required  a  more  intricate  design  and  would  be  difficult  to 
control  in  an  automated  sequence.  The  method  used  now  consists  of  water  spraying  inside  a  non¬ 
oxidizing  tank  which  acts  as  a  heat  exchanger.  The  water  is  sprayed  with  the  use  of  a  commercially 
available  nozzle,  with  a  90°  full  cone  spray  angle.  Initially  the  use  of  three  spray  nozzles  (Figure 
3.10)  was  employed  to  decrease  the  modes  of  failure  in  the  subsystem,  but  ultimately  the 
redundancy  was  eliminated  and  only  one  nozzle  is  used.  The  water  supplied  by  the  water 
connections  in  the  laboratory  fluctuates  from  15-20  gallons  per  minute  (gpm),  and  the  water 
pressure  is  70  psi.  The  spray  nozzle  is  design  to  supply  13  gpm  of  water  at  a  supply  pressure  of  70 
psi,  and  consists  of  a  0.31”  orifice  diameter. 


Figure  3.10:  Initial  conceptual  model  of  exhaust  section  with  spray  nozzle  diagram 
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The  analysis  for  the  exhaust  cooling  was  done  both  theoretically  and  computationally, 
employing  equations  of  heat  transfer,  and  a  multi-phase  model.  In  order  to  find  the  amount  of 
water  needed,  a  system  of  equations  is  used  and  includes  the  mass  fraction  of  liquid  water  and 
steam.  The  equations  are  solved  using  the  properties  of  water  and  steam,  as  well  as  the  initial 
temperature  of  800K  such  that- 


Qexhaust  Qwater  "h  Qsteam 


(jhCp (Th2  ^cS) eXhaust  ~  (piliqUid.Cp(T0  Tvater))water  +  (plcp(j'vapor  T-water)  +  rhtl)steam  (3.5) 

Qwater  ~  ^ liquid  T  T^-steam  (3.6) 

The  computational  model  used  the  result  of  the  mass  of  water  from  the  previous  equations 
at  a  temperature  of  293K,  and  pressure  of  70  psi.  The  total  mass  of  the  products  of  combustion  at 
800K  and  ambient  pressure  is  used.  A  transient,  multi-phase  turbulent  k-s  model  was  used  to  carry 
out  the  simulation  with  dimension  of  the  spray  nozzle  and  water  tank  as  specified  by  the 
manufacturer.  The  model  was  run  using  ANSYS  14.1  Academic  version  with  a  total  of  200000 
elements.  The  results  are  shown  in  Figure  3.11,  and  the  results  of  the  computational  temperature 
compared  to  the  theoretical  value  are  presented  in  Table  3.3.  The  minimum  water  flow  rate 
required  to  cool  the  exhaust  gas  is  7.5  gpm,  which  is  less  than  the  maximum  flow  provided  by  the 
spray  nozzle.  A  pump  is  used  in  order  to  drain  the  accumulated  water  out  of  the  tank  at  a  flow  rate 
higher  than  the  nozzle.  A  supplementary  dilution  line  was  also  added  to  the  exhaust  in  case  an 
unbumed  mixture  entered  the  exhaust. 


Table  3.3:  Comparison  between  theoretical  and  CFD  results 


Theoretical  exhaust  temperature 

CFD  mean  temperature 

333K 

340K 
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Figure  3.11:  Temperature  contours  in  exhaust  simulation 


3.4  Combustion  Chamber  Cooling 

As  temperature  rises  inside  the  combustor  during  combustion,  the  structural  integrity  of  the 
combustor  subsystem  is  compromised,  and  running  time  becomes  an  important  issue.  The  use  of 
various  cooling  methods  were  taken  into  consideration,  but  a  pre-existing  heat  exchanger  was 
rerouted  to  flow  coolant  through  the  combustor.  Figure  3.12  illustrates  the  conceptual  an  actual 
setup  of  the  heat  exchanger.  Figure  3.13  shows  additional  information  about  the  heat  exchanger 
which  can  supply  a  cooling  load  of  200kW.  In  order  to  address  the  temperature  rise  within  the 
combustor  during  the  combustion  process,  cooling  channels  were  manufactured  into  the 
combustor.  These  cooling  channels  are  9.5”  long,  have  a  0.5”  diameter  and  there  are  4  channels  in 
total.  These  channels  will  flow  dynalene  HC-10,  a  non-toxic  coolant,  at  a  rate  of  5.6  -  10  gpm. 
The  properties  of  this  coolant  are  obtained  from  the  manufacturer  [25].  The  temperature  at  different 
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critical  points  in  the  combustor  was  calculated  using  equations  for  transient  heat  transfer  on  a 
convective  surface  [26].  These  calculations  are  used  with  the  log  mean  temperature  difference 
method  (LMTD)  calculations  associated  with  the  coolant  flow  inside  the  combustor  channels. 


qs(P )  =  h[Tm  -  7X0,t)] 


(3.7) 


T(x,t)—T  i 
Too-Ti 


(3.8) 


A7)n  = 


ATe-ATj 


(3.9) 


These  equations  were  solved  for  different  flow  conditions  inside  the  combustor  and  inside 
the  cooling  channels  to  find  a  relationship  between  the  coolant  and  the  temperature  at  certain 
thicknesses  of  the  combustor  walls.  Equation  3.8  is  also  applied  to  the  quartz  windows  in  order  to 
calculate  the  difference  in  thermal  expansion  of  the  stainless  steel  and  quarts  interfaces  [27]. 


—  =  aAT  ;  —  =  2aAT  ;  —  =  3aAT  (3.10) 

L0  Aq  V0 

Comparing  the  thermal  expansion  differences  using  the  temperatures  from  Figure  3.14,  the 
adequate  gasket  sizing  was  chosen  to  be  1/16”  thickness. 


Figure  3.12:  CAD  model  of  heat  exchanger  and  dynalene  tank  (left)  and  physical  setup  (right) 
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Figure  3.13:  Additional  schematic  of  heat  exchanger  used  for  cooling 


Temperature  at  min  thickness  (no  cooling)  Temperature  at  min  thickness  (cooling)  Quartz 


Figure  3.14:  Temperature  vs  time  plots  for  critical  thickness  and  largest  quartz  window 


Using  the  thermal  expansion  of  stainless  steel  and  quartz  with  Figure  3.14,  the  determined 
allowable  temperature  is  125  °C.  The  allowable  run  time  is  doubled  when  the  combustor  is  being 
cooled  with  the  dynalene. 
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3.5  Control  System 


The  control  system  setup  was  designed  with  two  objectives  in  mind;  to  be  able  to  control 
all  the  instrumentation  manually  as  well  as  automatically  through  the  use  of  LabVIEW  with 
National  Instruments  PCI  cards,  and  to  be  efficient  and  simple  to  operate.  In  order  to  achieve  this, 
all  instrumentation  was  wired  into  a  single  space  consisting  of: 

Main  control  box  with  connector  blocks  attached  which  controls: 
o  AC  valves 
o  DC  valves 

o  Manual  emergency  stops 
o  Spark  plug  ignition 

DC  power  supplies  for  valves  and  instrumentation  (thermocouples,  transducers) 
LabVIEW  designated  CPU  for  running  automated  sequences 
Signal  generator  for  pilot  flame  ignition 

This  control  box  is  shown  in  Figure  3.15  and  the  corresponding  electrical  schematic  is 
shown  in  Figure  3.16.  There  are  two  separate  manual  emergency  stops  in  the  control  system;  one 
for  all  gas  controlling  valves,  and  one  for  water  controlling  valves.  The  control  box  also  consists 
of  a  double  pole-double  throw  (DPDT)  switch  which  serves  as  a  manual  override  in  case  of  a  CPU 
or  software  malfunction.  The  LabVIEW  program  consists  of  four  Data  Acquisition  assists  to 
accommodate  for  the  different  signal  each  type  of  instrumentation  requires.  Table  3.4  presents 
what  each  DAQ  assist  controls.  The  LabVIEW  interface  block  diagram  are  shown  in  Figure  3.17 
and  3.18  respectively. 


Table  3.4:  DAC  controlled  instrumentation 


DAQ  Assist 

Instrumentation  controlled/monitored 

1 

Pressure  transducers,  flowmeters,  and  thermocouples 

2 

Air/methane  line  valves 

3 

Air  and  methane  proportional  control  valves 

4 

Pilot  flame  valves 
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Figure  3.15:  Main  control  box 


NDTESi 

4.  ALL  DAO  RELAYS  (DR-01/DR-14)  ARE  POWERED  THROUGH  THE  DATA 
ACQUISITION  COMPUTER. 

3.  ALL  RELAYS  ARE  IN  THE  NO  POSITION.  ALL  RELAYS  ARE 

NON-LATCHING  AND  WILL  RETURN  TO  THE  ORIGINAL  POSITION  WVCN 
POWER  IS  REMOVED 
H>4  CHANNEL  DC  POWER  SUPPLY 

[j>PS-l  LABORATORY  WALL  POWER  OUTLET 


Figure  3.16:  Electrical  schematic  of  main  control  box 
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Figure  3.17:  LabVIEW  control  panel  interface 


Figure  3.18:  DAQ  assist  block  diagram 
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3.7  Test  Setup  Summary 


To  summarize,  the  system  can  be  described  using  Table  3.5.  The  subsystems  and  their 
function  are  explained.  For  the  present  work  however,  only  the  main  air,  exhaust,  and  PIV 
subsystems  are  considered. 


Table  3.5:  A  summarization  of  the  subsystems  that  make  up  the  high  turbulent  intensity  system 


Subsystem 

Function 

Combustor 

Provides  optical  access  in  order  to  employ 

flow  diagnostic  techniques  and  study  the 

previously  mentioned  research  topics 

Main  air  line 

Provides  cold/preheated  air  to  the  combustor 

for  premixed  flame  combustion 

Pilot  flame 

Ignition  source  for  the  air  and  methane 

reactants  inside  the  combustion  chamber 

Methane  line 

Provides  methane  to  the  combustor  for 

premixed  flame  combustion 

Exhaust 

An  exhaust  which  serves  as  an  escape  which 

cools  the  hot  combustion  products  to  a  safe 

and  acceptable  temperature 

Chamber  cooling 

Cools  the  combustor  in  order  to  prolong 

combustion  run  time  for  up  to  13  minutes 

Control  system 

Provides  manual  and  automatic  control  of 

system  instrumentation 

PIV  system 

Laser  diagnostics  system  used  in  the  present 

work  to  characterize  the  turbulence  scales  of 

air  at  compressible  conditions 

PLIF  system 

Laser  diagnostics  system  that  allows  the 

analysis  of  flame  front  characteristics. 
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Chapter  4:  Methodology 


The  present  work  focuses  on  the  study  of  turbulent  structures  in  the  flow  field  at 
compressible  air  flow.  Figure  4.1  shows  the  schematic  pertaining  to  the  setup  used  for  the  current 
work.  A  test  matrix  is  developed  to  test  the  three  different  blockage  ratios  and  volumetric  flow 
rates  as  seen  in  Table  4.1.  Before  the  tests  are  conducted,  the  PIV  laser  and  CCD  camera  are 
synchronized  and  calibrated  in  order  to  optimize  results.  All  data  acquisition  is  obtained  and  post- 
processed  using  the  ‘DynamicStudios’  software.  The  three  different  grids  used  are  described  in 
Table  4.2. 


NOTES 

1.  r  150  ANSI  Flange 

2.  5"  150  ANSI  Flange 

3.  1’  150  ANSI  Flange 

4.  3"  150  ANSI  Flange 


Figure  4.1:  Flow  schematic  of  the  present  study 
Table  4.1:  Test  matrix  for  fluid  flow  studies 


Blockage  Ratio  (%) 

Fluid  Temperature  (K) 

Test  Pressure  (PSIG) 

Flow  Rate  (scfm) 

Re 

54 

298 

0 

25 

5762 

54 

298 

0 

55 

12676 

54 

298 

0 

115 

26506 

61 

298 

0 

25 

5307 

61 

298 

0 

55 

11676 

61 

298 

0 

110 

24413 

67 

298 

0 

25 

4852 

67 

298 

0 

55 

10675 

67 

298 

0 

115 

22321 
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To  carry  out  the  PIV  experiments,  a  '15-1000  Dual  power'  527nm  Neodymium-doped 
Yttrium  Lithium  Fluoride  (Nd-YLF)  laser  with  a  pulse  energy  of  15mJ/pulse  was  used.  Laser  is 
redirected  at  45°  through  a  telescopic  lens  which  then  generates  the  light  sheet  used  to  illuminate 
the  alumina  particles  seeded  into  the  flow.  The  image  plane  was  captured  by  a  'Speed  sense  9070' 
CCD  camera  with  a  frame  rate  of  3.1  kHz.  The  image  plane  was  defined  by  the  area  the  camera 
captures  which  contains  the  laser  sheet;  this  area  is  a  64mm  x  40mm  space.  The  flow  field  was 
seeded  with  lpm  alumina  particles,  and  these  particles  are  seeded  far  upstream  to  where  the 
introduction  of  particles  does  not  affect  the  cross-flow  development.  The  camera  and  the  laser 
pulse  are  synchronized  and  controlled  using  a  BNC  575-8  time  box.  The  timer  box  uses  BNC 
cables  to  synchronize  the  laser  firing  pulse  and  the  camera  trigger  rate.  The  signals  of  both 
components  are  routed  through  the  timer  box  which  synchronizes  them  in  the  DynamicStudios' 
software.  The  image  is  calibrated  in  order  to  form  a  correlation  between  pixel  values,  and  physical 
dimensions.  The  image  is  calibrated  while  the  laser  is  firing  and  particles  are  present  in  the  flow 
to  optimize  the  camera  focus.  Ideally,  the  camera  will  be  focused  in  the  center  of  the  light  sheet  to 
capture  the  highest  laser  intensity  plane.  Raw  PIV  images  were  taken  at  a  rate  of  3  kHz  in  a  double 
frame  capture  mode,  and  a  total  number  of  375  images  were  used  for  analysis.  The  image  number 
is  obtained  from  the  time  step  defined  in  Chapter  2,  section  2.2.1.  A  raw  image  used  for  image 
calibration  is  shown  in  Figure  4.2  and  a  raw  image  detailing  the  image  characteristics  is  shown  in 
Figure  4.3. 


Figure  4.2:  Raw  PIV  image  showing  calibration 
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Figure  4.3:  Raw  image  detailing  the  different  components  in  the  PIV  capture 

The  main  air  line  is  controlled  by  three  different  valves.  A  quarter  turn  manual  control 
valve  is  used  to  allow  air  to  enter  the  main  airline  into  the  system  from  the  compressor  mentioned 
in  Chapter  3.  The  air  passes  through  a  120  V  (AVCO  E03/05  8P14)  solenoid  valve  which  has  a 
maximum  torque  of  35  Nm  at  0.5  amp  and  a  12  second  open/close  time.  After  that  the  flow  passes 
through  a  DC  ( EZ  88E22-C0EWS-D01 )  proportional  control  valve.  The  proportional  control  valve 
allows  a  0  -  10V  input  signal  to  control  how  much  the  valve  opens  to  allow  air.  The  air  passes 
through  an  analog  flowmeter  ( Hedland  variable  area  flowmeter )  which  has  flow  measurement 
range  of  25  -  250  scfm  at  a  pressure  range  of  40  -  120  psi  and  ±2%  accuracy.  Figure  4.4  shows 
the  solenoid  valve,  proportional  control  valve,  and  flowmeter.  Figure  4.5  shows  the  laser  system 
in  operation.  Figure  4.6  shows  an  overall  view  of  the  main  air  line  and  seeding  point  of  alumina 
particles. 


35 


Figure  4.4:  Solenoid  valve  (a),  proportional  control  valve  (b)  and  flowmeter  (c) 


Figure  4.5:  PIV  system  being  fired  into  the  combustor  for  flow  analysis 
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Figure  4.6:  Main  airline  left  and  mixing  chamber  with  seeding  point  indicated  (right) 


Table  4.2:  Grid  description 


Blockage  ratio  (BR) 

54% 

61% 

67% 

Effective  area  (in2) 

0.46 

0.39 

0.33 

Hole  area  (in2) 

0.0184 

0.0156 

0.0132 

Hole  diameter  (in2) 

0.153 

0.140 

0.129 

Hole  diameter  (mm) 

3.8 

3.5 

3.2 

37 


Chapter  5:  Result  and  Discussion 


As  previously  mentioned,  initial  system  validation  has  been  conducted  using  PIV  [28].  It 
was  concluded  that  the  flow  field  in  the  combustor  exhibited  the  same  behavior  as  seen  in  the 
LES  simulations.  The  velocity  profile  showing  agreement  between  the  two  methods,  is  shown  in 
Figure  5.1.  Additional  testing  was  done  in  the  present  work  to  further  validate  the  flow  field. 
Figure  5.2  shows  an  almost  identical  PIV  image  to  that  obtained  in  the  FES  model  for  Re  =  815. 


Figure  5.1 :  Average  u  velocity  profile  after  grid  comparison  between  CFD  and  PIV  results  [28] 


Figure  5.2:  comparison  of  PIV  images  and  CFD  results  for  RE  =  850 
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5.1  Particle  Image  Velocimetry  (PIV)  Testing 


The  results  of  the  experiments  carried  out  for  the  present  work  are  summarized  in  Table 
5.1.  Three  sets  of  images  were  generated  per  blockage  ratio  and  flow  condition,  equaling  a  total  of 
27  image  sets  for  analysis.  Of  those  image  sets,  further  analysis  narrowed  down  the  post¬ 
processing  to  9  image  sets,  with  375  images  per  analysis.  The  results  of  post-processing  are 
presented  in  the  next  sections  to  show  how  the  flow  structures  differ  per  grid  and  flow  rate.  The 
results  from  PIV  imaging  were  taken  at  a  line  before  the  step  and  one  after  the  step  (Figure  5.1)  to 
analyze  the  flow  structure  changes  caused  by  the  backward  facing  step.  From  the  images,  proper 
orthogonal  decomposition  (POD)  analysis,  kinetic  energy  measurement,  and  turbulent  intensity 
measurements  were  obtained. 


Table  5.1:  Summarization  of  flow  condition  results 


Blockage  Ratio  (%) 

Flow  Rate  (scfim) 

Re 

Proportional  Valve  Voltage 

67 

25 

4852 

3.5 

67 

55 

10675 

4.2 

67 

115 

22321 

5.7 

61 

25 

5307 

3.5 

61 

55 

11676 

4.2 

61 

115 

24413 

5.7 

54 

25 

5762 

3.5 

54 

55 

12676 

4.2 

54 

115 

26506 

5.7 
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Figure  5.3:  Area  of  interest  for  PIV  analysis  with  lines  indicating  location  of  data  extraction 


5.1.1  Proper  Orthogonal  Decomposition  Results 

Proper  orthogonal  decomposition  (POD)  is  a  post-processing  technique  offered  in 
‘DynamicStudios’  which  allows  PIV  data  to  be  decomposed  into  different  energy  modes.  First  a 
snapshot  is  taken  from  a  range  of  images  which  gives  a  graphical  representation  of  the  energy 
decomposition  in  the  flow.  The  mean  velocity  field  is  considered  the  zero'th  mode  of  the  POD. 
Subtracting  the  mean  from  all  snapshots,  the  rest  of  the  analysis  operates  on  the  fluctuating  parts 
of  the  velocity  components  (umn,  Vmn,  Wmn)  where  u,  v  &  w  denote  the  fluctuating  part  of  each 
velocity  component  [29].  The  fluctuation  components  are  then  resolved  into  a  multi-dimensional 
matrix  and  decomposed  into  the  different  energy  modes  in  the  flow.  Figure  5.4  shows  a  sample 
POD  snapshot,  Appendix  IV  has  a  complete  list  of  all  POD  snapshots.  This  snapshot  is  then 
projected  onto  the  PIV  images  to  isolate  the  flow  structure  according  to  the  chosen  energy  modes. 
The  POD  projections  for  the  three  blockage  ratios  are  presented  in  Figures  5.5-5.10. 
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Modal  Energy  Distribution 


Figure  5.4:  Sample  POD  snapshot  for  BR  67%  and  Re  =  22321 


When  carrying  a  POD  analysis,  the  POD  snapshot  is  crucial  to  the  validation  of  results. 
The  snapshot  should  exhibit  an  exponential  decrease  in  energy  as  mode  number  increases.  This 
indicates  the  concentration  of  energy  in  the  first  modes  which  are  associated  with  the  largest  scales 
in  the  fluid.  This  analysis  provides  a  quantification  of  the  integral  length  scale,  Lt,  defined  in 
Chapter  2.  The  lower  energy  modes  are  noise  and  therefore  not  of  interest. 
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Figure  5.5:  POD  projection  with  67%  BR  and  Re  =  22321 
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Figure  5.6:  POD  projection  with  67%  BR  and  Re  =  10765 


Figure  5.7:  POD  projection  with  61%  BR  and  Re  =  24413 
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Figure  5.8:  POD  projection  with  61%  BR  and  Re  =  1 1676 


Figure  5.9:  POD  projection  with  54%  BR  and  Re  =  26506 
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Figure  5.10:  POD  projection  with  54%  BR  and  Re  =  12676 

The  POD  analysis  shows  a  comparison  between  lower  and  higher  Reynolds  numbers,  and 
shows  the  structural  difference  in  the  energy  modes.  The  modes  represent  the  structures  which 
carry  the  highest  amount  of  energy.  The  largest  energy  concentration  regions  follows  a  similar 
pattern  for  all  cases.  In  the  lower  flow  cases  (Figures  5.6,  5.8  and  5.10)  the  flow  carries  the  highest 
energy  at  the  entrance  of  the  combustor  and  at  the  region  in  the  step.  Clear  recirculation  is  only 
seen  in  Figure  5.10  and  5.8  however,  and  in  Figure  5.6  there  is  a  shear  layer  which  cause  the  flow 
to  separate  after  the  step.  In  figure  5.6  the  appearance  of  jets  caused  by  the  perforated  plate  can  be 
seen,  and  these  jets  break  up  once  the  flow  passes  the  step.  In  the  higher  flow  cases  the  flow 
structures  are  dominated  by  the  entrance  regions  and  shear  layers  similar  to  that  shown  in  Figure 
5.6.  The  appearance  of  jets  is  apparent  on  Figures  5.5,  5.7  and  5.9  though  only  Figure  5.5  shows 
clearly  defined  flow  jets  which  carry  past  the  light  sheet  area.  Figure  5.7  also  shows  a  clear  jet  in 
the  center  of  the  flow  profile,  though  this  jet  breaks  up  before  the  end  of  the  light  sheet.  Though 
Figure  5.9  exhibits  favorable  breakup,  the  shear  layer  in  Figure  5.7  shows  a  reattachment  point 
toward  the  end  of  the  light  sheet. 


44 


5.1.2  Measurement  of  Kinetic  Energy 


The  kinetic  energy  distribution  is  compared  for  the  three  grids  at  the  lines  before  and  after 
the  step.  Equation  2.6  is  used  to  calculate  the  turbulent  kinetic  energy,  and  these  values  are  then 
normalized  when  plotted.  The  results  of  kinetic  energy  with  respect  to  y  -  position  of  a  particle  are 
shown  in  Figure  5.12  -  Figure  5.15.  The  normalized  kinetic  energy  will  serve  as  a  way  to  analyze 
how  the  energy  is  distributed  before  and  after  the  step,  as  well  as  to  show  how  this  energy  increases 
when  flow  rate  is  increased.  Figure  5.11  shows  a  general  visual  representation  of  kinetic  energy 
distribution  in  the  flow  field. 
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Figure  5.12:  Normalized  kinetic  energy  before  the  step  for  BR  =  54%,  61%  and  67%  at  Re  = 

12676,  11676  and  10675  respectively 
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Figure  5.13:  Normalized  kinetic  energy  before  step  for  67%  and  61%  BR  at  22321  and  24413 

respectively 
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Figure  5.14:  Normalized  kinetic  energy  after  the  step  for  BR  =  54%,  61%  and  67%  at  Re  = 

12676,  11676  and  10675  respectively 


Figure  5.15:  Normalized  kinetic  energy  after  the  step  for  BR  =  54%,  61%  and  67%  at  Re  = 

26506,  24413  and  22321  respectively 
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At  the  lower  flow  rate,  BR  =  67%  has  the  most  TKE  compared  to  the  other  two  grid 
geometries,  but  at  the  higher  flow  rates  the  TKE  favors  BR  =61%  both  before  and  after  the  step. 
Figure  5.14  shows  the  three  grid  TKE  values  follow  a  similar  trend  caused  by  the  grid  induced 
turbulence.  Much  like  in  the  POD  analysis,  the  BR  =  67%  shows  less  breakup  in  the  flow  as  seen 
by  the  drop  in  TKE  shown  in  Figure  5.15.  Again  at  higher  flow  rates  the  BR  which  shows 
promising  values  is  61%. 

5.1.3  Measurement  of  Turbulent  Intensity 

The  turbulent  intensity  is  an  important  parameter  in  determining  turbulence  in  the  flow, 
and  more  importantly  the  velocity  fluctuations  described  in  Section  2.4.  General  values  of 
turbulent  intensity  range  from  5  -  20%.  Turbulent  intensity  is  calculated  before  and  after  the  step 
and  compared  for  the  different  grid  geometries  and  flow  rates.  Equation  2.5  is  used  to  calculate 
the  value  of  turbulent  intensity.  The  results  for  turbulent  intensity  are  presented  in  Figure  5.16  — 
Figure  5.19. 


Figure  5.16:  Turbulent  intensity  before  step  with  BR  =  54%,  61%  and  67%  and  Re  =  12676, 

1 1676  and  10675  respectively 
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Figure  5.17:  Turbulent  intensity  before  step  for  BR  =  67%  and  61%  with  Re  =  24413  and  22321 

respectively 
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Figure  5.18:  Turbulent  intensity  after  the  step  with  BR  =  54%,  61%  and  67%  and  Re  =  12676, 

1 1676  and  10675  respectively 


49 


0  5  10  15  20  25  30  35 

y-position  (mm) 

Figure  5.19:  Turbulent  intensity  after  the  step  with  BR  =  54%,  61%  and  67%  and  Re  =  26506, 

24413  and  22321  respectively 

Again  BR  =  67%  exhibits  greater  turbulent  intensity  at  lower  flow  rate  values,  but  drops  in 
the  case  of  maximum  flow.  Figure  5.16  and  Figure  5.18  show  the  BR  =  67%  turbulence  intensity 
is  greater  than  the  other  grids  which  validates  this  assumption.  Figure  5.19  shows  how  the  turbulent 
intensity  in  BR  =  67%  drops  below  the  other  grids  and  BR  =  61%  is  the  dominant  turbulence 
generator;  this  assumption  is  also  validated  in  Figure  5.17.  It  should  be  noted  that  in  Figure  5.17 
the  turbulent  intensity  follows  a  parabolic  profile  which  has  a  peak  at  the  bottom  and  top  surfaces 
of  the  combustor.  This  behavior  can  be  attributed  to  the  formation  of  jets  seen  in  the  POD  analysis 
for  both  BR  =  67%  and  61%.  This  parabolic  profile  is  not  apparent  in  Figures  5.16  and  5.18, 
correlating  to  the  formation  of  the  jets  at  higher  flow  rates  shown  in  the  POD  analysis.  Figure  5.19 
illustrates  a  more  predictable  trend  for  turbulent  intensity  in  the  grids  meaning  the  jets  begin 
breaking  up,  causing  fluctuations  in  the  flow  field. 


50 


5.1.4  Comparison  of  Grid  Thickness 


From  the  previous  analyses,  it  was  determined  that  the  formation  of  jets  affects  the 
behaviors  of  the  flow  and  causes  the  grid  to  act  as  an  orifice  rather  than  a  turbulent  generator.  In 
order  to  address  this  issue,  the  BR  =  67%  grid  was  duplicated  and  its  thickness  was  halved.  This 
BR  was  chosen  due  to  being  the  most  affected  by  the  jet  formation  in  the  previous  sections.  Figure 
5.20  shows  a  POD  projection  for  the  thinner  grid,  Figure  5.21  shows  a  comparison  of  TKE,  and 
Figure  5.22  shows  a  comparison  of  turbulent  intensity.  The  purpose  of  changing  the  grid  thickness 
was  to  determine  whether  significant  improvements  are  made  to  the  flow  structure  in  order  to 
warrant  a  more  detailed  analysis  of  the  grid  geometry. 


Figure  5.20:  POD  projection  of  thin  grid  BR  =  67%  and  Re  =  22321 

Comparing  Figure  5.20  to  Figure  5.5  shows  a  dramatic  change  in  flow  structure.  The  jets 
begin  breaking  up  much  sooner  in  the  thin  grid,  and  there  is  some  recirculation  apparent  under  the 
shear  layer  after  the  step.  The  recirculation  in  the  steps  increases  the  value  of  the  integral  length 
scale  which  benefits  the  scope  of  the  present  work.  Previously  the  length  was  determined  by  the 
size  of  the  jets,  which  were  assumed  to  be  the  size  of  the  perforations  in  the  grid. 
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Figure  5 .2 1 :  Normalized  kinetic  energy  for  regular  and  thin  grid  at  BR  =  67%  and  Re  =  2232 1 


Figure  5.22:  Turbulent  intensity  comparison  between  regular  and  thin  67%  BR  at  Re  =  22321 

Figure  5.21  and  Figure  5.22  show  the  values  of  TKE  and  turbulent  intensity  do  not  vary 
dramatically  with  the  change  in  grid  thickness.  Figure  5.21  shows  a  slight  increase  in  TKE  for  the 
thinner  grid  which  can  be  attributed  to  the  earlier  breakup  of  the  jets.  Figure  5.22  however  shows 
the  turbulent  intensity  remains  unchanged  for  the  two  grids.  The  results  of  changing  the  grid 
thickness  suggest  a  further  investigation  of  grid  geometry  should  be  conducted  to  find  an  optimal 
grid  geometry. 
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Chapter  6:  Conclusion  and  Future  Work 

6.1  Summary 

Flow  mapping  of  a  high  turbulent  intensity  combustor  has  been  completed.  An  analysis  of 
high  Reynolds  number  fluid  flow  over  a  backward  facing  step  has  been  performed  using  PIV 
imaging.  Results  of  the  analysis  indicate  the  grid  geometry  has  a  great  effect  on  flow  structure  as 
evidenced  by  the  results;  the  conclusion  being  that  a  thinner  grid  with  blockage  ratio  =  61% 
provides  better  flow  breakup.  These  results  confirm  the  presence  of  high  turbulence  in  the  system, 
and  the  varying  length  scales  associated  with  different  geometries.  The  three  main  parameters  were 
resolved  for  different  turbulent  flow  rates,  including  compressible  flow.  These  three  main 
parameters  serve  to  further  the  understanding  of  flow  interactions  and  how  these  affect  the 
turbulent  parameters  required  to  achieve  a  thickened  flame. 

POD  analysis:  provides  a  method  of  isolating  the  energy  modes  in  a  flow  to 
determine  structural  parameters  (Lt) 

Kinetic  Energy:  provides  additional  flow  mapping  and  denotes  the  profiles  in 
which  the  largest  energy  is  carried  by  the  flow 

Turbulent  Intensity:  Detennines  the  ratio  of  turbulence  affecting  the  flow  field  in 
order  to  quantify  velocity  fluctuations  in  a  flow  (u’) 

These  parameters,  along  with  combustion  parameters,  will  determine  the  flame  regime 
studies  for  the  high  turbulent  intensity  system. 
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6.2  Future  Work 


Further  grid  studies  can  be  performed  in  order  to  find  an  optimization  point  for  grid 
blockage  ration  and  thickness.  Furthermore  the  combustion  system  is  ready  to  begin  testing  and 
analyzing  flame  behavior.  The  further  considerations  of  the  project  can  be  described  as  follows: 

To  further  the  development  of  a  backward  facing  step  stabilized  flame  system  by 
conducting  initial  flame  studies  using  PLIF. 

To  study  flow  and  flame  characteristics  using  simultaneous  Time  Resolved 
Particle  Image  Velocimetry  (TR-PIV)  and  OH/CH  Planer  Laser  Induced 
Fluorescence  (PLIF)  imaging. 

To  define  the  effects  of  bulk  velocity  parameters,  turbulence  parameters  and 
recirculation  zone  characteristics  on  flame  location,  shape,  fluctuation,  reaction 
zone  and  local  displacement  speed. 

The  findings  presented  in  this  work  will  be  used  with  the  future  flame  studies  in  order  to 
determine  and  characterize  a  thickened  flame.  The  experimental  findings  of  a  flame  in  the  thin 
reaction  zone  regime  will  provide  useful  data  to  develop  numerical  and  computational  models 
which  can  be  applied  to  the  design  of  high  turbulent  combustion  systems. 
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Appendix  I 

Specification  of  pressure  transducers  and  thermocouples. 


LORD 


STELLAR  TECHNOLOGY 


I  t- 02 


?J7  Ccmm**c#  0**** 

Hew  VerK  1422*  LSA 
25CL19CO  •  r«x:Tid230  1«0B 


TEMPERATURE  CALIBRATION  RECORD 

CUSTOMER 

SERIES*>ART  NUMBER 

SERIAL  NUMBER 

SALES  ORDER 

University  Of  To 

*»e  1?  El  V*»o 

PT17OO-150G-171 

240144$ 

14*0026 

CUSTOMER  PO 

CUSTOMER  PART  NUMBER 

MAXiCAL  VOLTAGE 

2024047243 

N/A 

10/12  V<k 

TEMP  MANGE  fl'l 

OUTPUT  RANGE 

CONDITION 

TO  300 

0-5  V 

NEW 

AU  DATA  IN  ’F 


STEP 

TEMP/TI 

Output.  V 

litul.  V 

Error,  V 

Error.  %  FSO 

0.01 

0.003 

0.001 

0.000 

0.001 

0.02 

12.  *>% 

CfjA 

0.640 

0 . 626 

0.015 

0.301 

2S  OS 

75.000 

1.275 

1.250 

0.025 

0.496 

37.5% 

112.500 

1.00$ 

i.873 

0.031 

0.615 

50.0% 

150. 000 

2.333 

2.300 

0.033 

0.654 

«•« 

1B7.50O 

3.156 

3.125 

0.031 

0.615 

r  vi.ot 

225-000. 

3.775 

3  760 

0.025 

0.497 

«1.5% 

.’62.5U0 

4.390 

4  .  175 

0.015 

0.301 

1  100.0% 

JUO.DQO 

5 .  c  a  i 

5.000 

0.001 

0.326 

ERROR  DISTRIBUTION 


SUMMARY 

CONNECTIONS 

RISC 

VALIDATION  POtNT  149  <TF  2 . 51 6 

PIN 

DESCRIPTOR 

CONNECTOR 

STATIC  ACCURACY  l\  FSO)  0.654 

A 

ESC- 

202031 

D 

EXC- 

CM  NEC  NO 

B 

PSIG* 

41912.42100  I 

Q 

PSIG- 

QA  STAMP 

t 

rsiGt- 

F 

rsiG- 

(  48  J 

V2*y 

9/30/2C14 

sa%mm  jam* 
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LORD  STTT  SfELLAR  technology 

137  Commarc*  Oma 

Anrtttti  Mea  von  1-U2IUSA 

T*T  7T4.7 W  T wo  .  F.I714.M01W* 

TEMPERATURE  CALIBRATION  RECORD 

CUSTOMER 

SERlESiPART  NUMBER 

SERIAL  NLTMBER 

vursoROHR 

University  of  7r*b«  6  Cl  Ease 

prl700-150G-171 

1401447 

[14-0026 

CUSTOMER  PO 

CUSTOMER  PART  NlMBEH 

MAXCAl  VOLTAGE 

2014047243 

N/A 

18/12  Vdc 

TEMP  RANGE  <T) 

OUTPUT  RANGE 

CONDITION 

TO 

V 

NEW 

ALL  DATA  <H  T 


STEP 

tempcfl^ 

OutpcA  V 

IdoM,  V 

tm if.  V 

error  %  ESO 

0.000 

-0*003 

0.000 

-0 .003 

-0.068 

12.5% 

37.500 

0.635 

0 . 62  5 

0.010 

0.207 

25.0% 

75.000 

1.2TO 

1.250 

0.020 

0.403 

37.5% 

112. 50C 

1.901 

1.875 

0.026 

0.521 

50.0% 

150.000 

2.S2E 

2.50C 

0.028 

0.560 

€2.5% 

107. son 

3.151 

3.125 

0.026 

0.521 

T5.0» 

22S.ODO 

3.770 

3.750 

0.020 

C.403 

•9.  At 

262. S00 

4.385 

4.375 

o.  gio 

0.207 

lou.u% 

lUO.OOO 

4.997 

5.000 

-0.0C3 

-0. 06? 

ERROR  DISTRIBUTION 


*  fc— A  %  TVS 


SUMMARY 

COWSYcnows 

MI9C 

VALUATION  PONT  1507TF  2.528 

PIN 

DESCRIPTION 

CONNECTOR 

STATIC  ACCURACY  |*FS0|  0.560 

A 

LXO 

202031 

D 

EXC- 

CA1  REC  NO 

B 

P3IG+ 

41912. 1«*B9 

C 

PS1G- 

QA  STAMP 

E 

TS1G+ 

r 

TS1G- 

9730/201 4 
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■  .  237  Com»«rc«  Oflv* 

LORD  S7X  sifiub  TFCHNotoGv  f-| -j\  rssssxsr.'sssm.m. 


CALIBRATION  RECORD 

CUSTOMER 

University  of  Texas  9  El  Paso 

SERIES.PART  NUMBER 

FT170D-1SQC-171 

1401446  14-0026 

CUSTOMER  PO 

2014041243 

CUSTOMB»  PART  NUMBER 

M/A 

PRESSURE  RANGE 

0-150  MIG 

COMP  TEMP  RANGE  <T) 

TO 

OPER.  TEMP  RANGE  (T) 

TO 

GAGE  RESISTANCE 

5000 

MAXiCAL  VOLTAGE 

VOLTS 

PROOE  i  BURST  PRESSURE 

225/30C  PSI 

CONDITION 

NEW 

ALL  DA r A  IN  V  TEST  TEMPERATURE:  AMBIENT 
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LORD  STC  Sll  I  I  \U  II  (.  IINOKK.V 


fT  -o-l 


737Co»nmcicc  OOve 
MWII  Hrti  Yotk  14 22V  USA 
T«l  nt  23D.120D  • 


CALIBRATION  RECORD 

CUSTOMER 

- KRTCSMARfNUM'bEN 

SERIAL  NUMBER 

SALES  ORDER 

University  of  Texas  0  El  Paso 

P71700-150G-171 

140144? 

14-0026 

■CUSTOMER  P0 

CUSTOMER  PART  NUMBER 

PRESSURE  RANGE 

2014047243 

M/A 

0-150  PS1G 

COMP  TEMP  RANGE  |*F| 

GAGE  RESISTANCE 

PROOF !  BURST  PRESSURE 

TO 

5000 

225/ 3D 0 

PSI 

OPCR  TIME  RANGE  |*T| 

MAXiCAL  VOLTAGE 

CONDITION 

TO  230 

IB/12  VOLTS 

?<£M 

ALL  DATA  IN  V  TEST  TEMPERATURE:  AMBIENT 


INCREASING 

DECREASING 

STEP 

OUTPUT 

BEST 

DIFFERENCE 

OUTPUT 

BEST 

DIFFERENCE 

HYSTERESIS 

0% 

-0.008 

-0.006 

-0.002 

-0.006 

-0.006 

0.000 

0.002 

20% 

D.9B9 

0.990 

-0.001 

0.992 

0.990 

0.002 

0.003 

1.9*1 

1.901 

0.000 

1.989 

1.987 

0.002 

0.002 

on 

A  .  8=54 

2.904 

0.000 

2.995 

2.984 

0.001 

0.001 

<KTA 

U.sBiL 

3.360 

0.000 

3.  3bi 

3.980 

C.OjI 

0.CQ.  , 

100% 

4.476 

4. 971 

-0.001 

4.976 

4 . 917 

-0.001 

o.ooc  i 

ISO 

4.994 

J 

SHUNT 

LOW  TEMPERATURES  HIGH  TEMPERATURE  2SITF 


STEP 

NORMALIZED 

STEP 

NORMALI2EO 

0% 

0.036 

OH 

0.063 

IOC f*, 

4 . 9?9 

100% 

4.992 

FSO 

4.943 

FSO 

4.929 

REPEATABILITY  fl  AMBIENT 


STEP 

OUTPUT  1 

OUTPUT  2 

0% 

-0.010 

-0.009 

100% 

4.974 

4.9"?i> 

FSO 

4-984 

4.934  | 
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Pulse  width  [ns]  Power  [W] 


Appendix  II 

Specification  of  PLIF  laser  and  Image  Intensifire. 

Laser  power  sheet 


Power  vs.  current 


Pulse  width  vs.  current 
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Pulse  energy  vs.  current 


0“ 

75- 

f=3KHz  50ns  |/%| 

f=5kHz  50ns 
f=10kHz50gs  |/\| 

7- 

6,5- 

6- 

- 

55- 

5- 

45- 

4- 

3,5- 

3- 

/ 

2  5- 

2- 

1  5- 

1- 

05- 

0- 

3 

5  36  37  3 

8  39  4 

0  41  42  43  4 

4  45  4 

6  4 

7  4 

8  4 

9  5 

0  5 

5 

2  5 

3  54 

Current  [A] 


I  frequency 

35  A 

- 40 - 

- 40 - 

50  A 

- 54S - 

0.00  w 

0,37  W 

7,53  W 

18.10  W 

23,47  W 

f=5kHz  50gs 

0,00  w 

0,18  W 

5,30  W 

17^95  W 

26,68  W 

fr  10kHz  50ns 

0,00  w 

0,00  W 

2,35  W 

1 1 .38  W 

20,52  W 

Frequency 

35  A 

40 - 

40 - 

50  A 

54  A 

f=3kHz  50ns 

9,4  ns 

7,0  ns 

6,0  ns 

5,7  ns 

f=5kHz  50ns 

8,7  ns 

7,0  ns 

6,5  ns 

f=10kHz  50ns 

12,7  ns 

9,1  ns 

7,8  ns 

Frequency 

50 - 

40  A 

40 

50  A 

54  A 

f=3kHz  50ns 

0,01  mJ 

0,18  mJ 

2,56  mJ 

6,32  mJ 

7,82  mJ 

f=5kHz  50ns 

0,01  mJ 

0,01  mJ 

1,08  mJ 

3,38  mJ 

5,34  mJ 

f=10kHz  50ns 

0,01  mJ 

0,01  mJ 

0,42  mJ 

1,19  mJ 

2,05  mJ 
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•Image  intensifier 


Parameter 

C10880-01 C 

C10880-03F 

Unit 

Luminous  sensitivity  (Typ.) 

230 

\ifij  Im 

Radiant  sensitivity  (Typ.)  15 

53 

mA/W 

Quantum  efficiency  (Typ.)  1 1 

15 

% 

Photocathode  window  material 

Synthetic  quartz 

- 

Photocathode  material 

Multialkali 

- 

Spectral  response  range 

185  to  900 

nm 

Wavelength  of  maximum  response 

430 

nm 

Photocathode  diameter 

25  4) 

mm 

Input  lens  mount 

C-mount 

F-mount 

Phosphor  screen  window  material 

G!  G33 

- 

Phosphor  screen  material 

P-46 

- 

Luminous  gain  (Typ.) 

Ixio5 

(Im/m2)/lx 

Radiant  emittance  gain  (Typ.)  1! 

4.2x103 

(W/m2)/(W/m2) 

EBI  (lumen)  (Typ.)  2) 

2x1 0'9 

Imfcm2 

EBI  (radiant)  (Typ.)  11  2) 

2x10'u 

W/cm2 

Limiting  resolution  (Typ.) 

38 

Lp/mm 

Image  magnification 

0.67 

- 

Maximum  input  light 
intensity  (Typ.) 

3> 

s.oxio'3 

lx 

1)  3) 

8.0x10'’° 

W/cm2 
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Appendix  III 

Sample  calculations 

In  order  to  prevent  choking,  we  analyze  the  required  cross  sections  needed  for  the  piping 
section.  Since  the  minimum  operating  pressure  is  3  bar,  then  the  upstream  pressure  must  be 
calculated  to  find  what  the  minimum  value  of  pressure  along  with  flow  rate  must  be  after  the 
compressor.  The  pressure  before  the  chamber  can  be  found  by  treating  it  as  an  isentropic, 
compressible  flow  in  a  nozzle.  From  literature: 

Air  at  M  =.  5 


Pt 

—  =  .8483  ->  P0  =  .8483/300000  =  353648 Pa 
°o 

This  can  be  used  as  the  final  pressure  after  the  heating  section.  Next  we  calculate  the  pressure 
drop  after  the  heating  section  (assumed  to  be  3m  max) 

Initially  at  the  compressor: 


/  m  \ 

Vmax  =  190  cfm  .  08976  —  ,  Pmax  =  140 psig  (9.853  bar), 
\s 


T  —  300 K  (slight  heating),  Ma  <  1 
Vmax  .08976 

v  =  =  _ =  457.873m/s 


A  .000196 


Since  the  velocity  along  the  pipe  is  higher  than  Ma=l,  then  flow  is  choked. 

To  prevent  this,  then  the  flow  after  the  heating  section  must  be  considered  since  it  will  have  a 
higher  velocity  than  that  at  the  inlet. 

After  heating: 


p  =  3.925  kg/m3  cp  =  1032  J/kgK 

Pr  =  .67  K=.0415w/mK 


p  =  2.707x10  5  kg/m*s  A  =  .000196m2 

g  =  9.81m/s2  p=  1/396.5 


m  .344715 

=  —  = - =  448.089m/s 

pA  3. 925(. 000196)  7 
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448.089 


Ma  =  —  =  - =  .99 

c  Vl-4(287)(500) 

We  want  Ma  <1  (.8  max); 


d  = 


4m 


InpMayJkRT 


.7 


A  diameter  of  0.7”  or  greater  will  prevent  choking. 


Fanno  line  flow  equations: 


T0 


=  T0 


* 


(fc+i)/ 

P0  _  p0  _  1/2  +  (k-  1  )Ma2\  Mk-i) 
P0*  p0*  Ma  \  k  + 1  ) 


f_=Po_  =  J_  f  k  +  1 

P*  ~  p0*  ~  Ma  V2  +  (k  -  1  )Ma2 


]^=Po_=  (  k  +  1 
T*  ~  p0*  ~\2  +  (k-  1  )Ma2 


V  _p*  _  If  k  +  1 
V*  ~  ~p  ~  ~Ma  V^+T^TjMa7 


fL*  _  1  -  Ma2  k  +  1  /  (k  +  1  )Ma2  \ 

D  kMa2  +  2/e  ^  \2  +  (/c  —  l)Ma2/ 
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Appendix  IV 

POD  snapshots 


Modal  Energy  Distribution 

Fraction  of  Total  Energy 


Fraction  of  Total  Energy 


POD  snapshot  for  BR  67%  Re  =  22321 

Modal  Energy  Distribution 


Fraction  of  Total  Energy 


POD  snapshot  for  BR  67%  Re  =  10675 


67 


Modal  Energy  Distribution 


Fraction  of  Total  Energy 


POD  snapshot  for  BR  61%  Re  =  24413 


Modal  Energy  Distribution 

Fraction  of  Total  Biergy 


POD  Mode  Number 


Fraction  of  Total  Energy 


POD  snapshot  for  BR  61%Re=  11676 
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Modal  Energy  Distribution 


Fraction  of  Total  &ergy 


POD  Mode  Number 


POD  snapshot  for  BR  54%  Re  =  26506 

Modal  Energy  Distribution 


Fraction  of  Total  &ergy 


POD  snapshot  for  BR  54%  Re  =  12676 


Fraction  of  Total  Energy 


Fraction  of  Total  Energy 
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Glossary 

(In  order  of  appearance) 


Symbol 

Re 

U 

L 

v 

u 

u' 

Urms 

i 

TKE 

n 

£ 

U\ 

Ilf] 

v 

>1 

U\ 

Rer 

Da 

Ka 

Lt 

St 

Sl 

M 

Q 

m 

cp 

As 

h 

T 

a 

k 


Description 
Reynolds  number 
Characterisrtic  velocity 
Characteristic  length 
Viscosity 

Instantaneous  velocity 
Fluctuation  of  velocity 
Root  mean  square  value  of  velocity 
Turbulent  intensity 
Turbulent  kinetic  energy 
Kolmogorov  length  scale 
Energy  dissipation 
Kolmogorov  time  scale 
Kolmogorov  velocity 
Fluid  viscosity 
Kolmogorove  length  scale 
Kolmogorove  time  scale 
Turbulent  Reynolds  number 
Damkohler  number 
Karlovivitz  number 
Integral  length  scale 
Flame  thickness 
Faminar  flame  speed 
Macn  number 
Heat  flow  per  unit  mass 
Mass  flow  rate 
Specific  heat 
Surface  area 

Convection  heat  transfer  coefficient 
Temperature 

Thermal  expansion  coefficient 
Thermal  conductivity  constant 
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